function dwdt=ifun2(t,w)%b e不等于0
q=1.6e-2;
e=1;
b=0;

m=0.02;%设置条件
dwdt=[w(2);
    q.*b./m.*w(4);
    w(4);
    (q.*e./m)-(q.*b./m).*w(2);
    w(6);
    0;];%建立微分方程
end